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We present results of a high precision Monte Carlo simulation of dynamically triangulated random surfaces (up 
to w 34,000 triangles) coupled to one scalar field (c = 1). The mean square extent has been measured for different 
actions to test the universality of the leading term as a function of the size of the surfaces. Furthermore, the 
integrated 2-point correlation function for vertex operators is compared with conformal field theory and matrix 
model predictions. 
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1. INTRODUCTION 

This talk reports about the results of several 
projects, where - apart from the authors - also M. 
Agishtein, R. Ben-Av, I. Klebanov, A. A. Migdal 
and S. Solomon were involved. Additional de- 
tails and even more precise data will be published 
soon. 

c = 1 matter coupled to 2-dimensional quan- 
tum gravity seems to be a critical model in the 
sense that most approaches to the theory of mat- 
ter coupled to gravity in 2 dimensions predict 
qualitatively different behaviour for c < 1 and c > 

1. Furthermore, there are two different theoreti- 
cal approaches (conformal field theory (l]||,||,[|,|| 
and matrix models |^,^J|,^,|l^,^l) ) which give 
quite detailed predictions for this model. 

The main aim of our high precision Monte 
Carlo simulations was to confirm those predic- 
tions where conformal field theory and matrix 
models agree, and to decide or clarify the situ- 
ation where the two approaches differ. 

2. DEFINITION OF THE MODEL AND 
THEORETICAL PREDICTIONS 

Let TJv denote the set of abstract triangula- 
tions of the 2-dimensional sphere with N trian- 
gles. This is equivalent to the set of planar, regu- 
lar graphs of degree 3 with N vertices and without 
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non-trivial 2-point subgraphs. Each graph can be 
characterized by its adjacency matrix: Cy = 1 if 
i is a neighbour of j and CVj = otherwise. The 
partition function of the models we simulated is 

Z N = J dX ^ ■■■ dX N (!) 

T~N / \ 

ex p( -lH c ij E (Xi,Xj) 

Xi is a real, scalar field attached to the faces of 
the triangulation. We mainly considered 



E 2 (X,Y) 



(2) 



the discrete analog of the squared derivative, but 
for a check of universal properties of some of the 
quantities we also present data for 



E 1 (X,Y) = \X-Y\ 



(3) 



Both actions are expected to describe c = 1 mat- 
ter coupled to quantum gravity. 

The 2-point function for the tachyon with mo- 
mentum p is 
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Gn( P ) = ^ 



N 2 ^ 



N ^ 



(4) 
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From this one obtains the moments of X by 
derivatives, e.g. the mean square extent 



(X 2 )n 




p=0 



(5) 



(G) 



For these quantities exist theoretical predictions 
from conformal field theory (l],|^,|],[||5| and from 
matrix models fl6]]7|J||^, |io||ll| ] ) , which can be sum- 
marized in the following formulas: 



G(p) 



1. Both approaches agree that 
z 3 



sinh 



(7) 



2. Furthermore, both approaches predict the 
general dependence of the mean square extent on 
the number of triangles: 

(X 2 ) N oc (lniV) 2 . (8) 

This equation, which has first been conjectured 
in (TJ|, implies that Gn(p) becomes a function of 
(plniV) only. 

3. In order to fix the coefficient in (||) one has 
to normalize the field X. This is usually done 
by using the asymptotic behaviour of the 2-point 
funcion in flat space: 



{X a X a ,) 
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For the Gaussian action (| 
known: 

a'{E 2 ) = 1/2tt « 0.159 



this normalization is 



(10) 



To determine this factor for action E\ of eq. (||) 
we performed a Monte Carlo simulation on a flat 
2-dimensional lattice (a 512 x 512 square lattice 
with periodic boundary conditions) and found 



a'(£i) « 0.134(2) 



(11) 



The conformal weight of the tachyon operator in 
flat space is 



A (P) = > 2 



(12) 



Using the KPZ formula [p[ to relate the weight 
of a conformal field in flat space with the one 
coupled to 2-dimensional gravity, 



A(p) 



24 A (p) - VT^5 



(13) 



one can apply formal scaling arguments to obtain 
G{p) as a function of the area A H|,|@: 



G A (p) oc A l - 2A & 



(14) 



For c = 1 this expression is non-analytic at p = 
0, and we get a wrong minus sign if we try to 
obtain the mean square extent directly from this 
formula, according to eqj^. For p large however, 
equation (Q) , which has been derived in conformal 
field theory by introducing a cut-off ||, agrees 
with (|l4|), provided that the mean square extent 
is given by 



(X 2 ) = ^(InAO 2 



(15) 



On the other hand, the mean square extent has 
been calculated in matrix models 0,^1) , yielding 



(X 2 ) = J (In TV) 2 . 

This calculation uses the action 



E a (X, Y) 



\X-Y\ 



(16) 



(17) 



Matrix models, however, do not only generate 
graphs corresponding to triangulations but also 
graphs with nontrivial 2-point subgraphs. This 
might imply that the "mass" has to be 

renormalized before one can compare the results. 

3. DATA AND RESULTS 

We made simulations for N =140, 420, 1260, 
3780, 11340 and 34020 triangles. The update 
of the scalar field was speeded up considerably 
by the use of the VMR ( "valleys-to- mountains 
reflections) cluster algorithm |TJJ]. For the up- 
date of the triangulation we used the triangle flip 

HJOIIDJUIl, which had t0 be su pp lemented 

by a local update procedure for the scalar field 
(we used both heat bath and Metropolis). Care- 
ful analysis of the autocorrelations revealed an 
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autocorrelation time of w 300 local sweeps (with 
a few clusters after each sweep) for the largest 
system, for which we performed about 200.000 
sweeps. 

No errorbars in the figures indicates that the 
errors are small compared to the symbols used to 
mark the data points. We should emphasize that 
the data for the simulations of action E\ are less 
accurate since we have not finished our runs. 
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Figure 1. Mean square extent for E 2 = \{X — Y) 2 



Fig. |l| shows the mean square extent for E2 as 
a function of N. A least square fit to 

< X 2 > N = a"(kiN) 2 +PQnN) +7 (18) 

yields a" = 0.080(4). The value for N = 8 can 
be calculated analytically and has been added for 
curiosity. 

Fig. ^ shows the corresponding data for action 
Ei, with a coefficient a" = 0.11(1) from a least 
square fit. While both actions quite nicely exhibit 
the (In ./V) 2 behavior, the coefficients are different, 
even if we take into account the normalization of 
the X-fields. This indicates that the coefficient 
might not be universal. The disagreement with 
the matrix model predictions ([l6]) can be due to 
the subtraction of graphs with non-trivial 2-point 
subgraphs. 
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Figure 2. Mean square extent for E\ = \X — Y\ 



Figs. H and [| show G(p) for the actions E 2 
(iV=34020 triangles) and E x (iV=3780 triangles) 
respectively, as well as a one-parameter fit with 
respect to (X 2 ) using the form (^). The univer- 
sality of the scaling function is nicely confirmed. 
For smaller lattices the curves, as a function of 
plnN, almost are on top of each other, which 
again implies the universality of eq. |[ Small de- 
viations from the theoretical prediction are not 
visible in the figures, but can be seen by compar- 
ing the numbers. They are probably due to finite 
size effects. From the values of (X 2 ) obtained 
by fitting G(p), the coefficient of (IniV) 2 is, con- 
sistently with the direct measurement presented 
before, 1.42 times larger for the action E\ than 
for E 2 - Whether we take into account the flat 
lattice normalization or not, this result indicates 
non- universality of a" . 

4. CONCLUSIONS 

The Monte Carlo simulations we performed for 
the mean square extent and the tachyonic 2-point 
function G(p) of c = 1 matter coupled to 2- 
dimensonal quantum gravity confirmed the uni- 
versal scaling behavior for G(p). Furthermore, 
they confirmed the (X 2 ) oc In 2 (TV) law and its 
universality. The coefficient in front of the In 2 
however seems not to be universal. The contra- 
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diction with matrix model predictions for this co- 
efficient might be due to a renormalization of the 
parameters in the matrix model. The discrep- 
ancy between two c = 1 scalar fields, which have 
the same conformal weight in flat space, can be 
a more serious problem. Presently, we perform 
several numerical checks to clarify these issues. 




Figure 4. G(p) for E 1 = \X- Y\ 
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